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. We study the fluctuations of the total topological charge of a scalar wave propagating 

1 I in a hollow conducting wave guide fllled with scatterers inside. We investigate the 

■ dependence of the screening on the scattering mean free path and on the presence 
' of boundaries. Near the cut-off frequencies of the wave guide, screening is strongly 

^ ' suppressed near the boundaries. The resulting huge fluctuations of the total topological 

^ ' charge are very sensitive to the disorder. 

•S! 

o : 

^ . 1. Introduction 
■ 

S ■ A complex scalar random wave field is given by: ^^(r, t) = A(r, t) exp(i$(r, t)) = 

\ x{r,t) + ir]{f,t). When both its real part x imaginary part 77 cancel, the 

I amplitude A also cancels but the phase is left undefined: it is a phase singularity. 

^ ■ In space, these singularities constitute nodal lines located at the intersection of 
the two surfaces defined by xi^^t) = and r](f^t) = 0. On a flat surface inside 

^ ■ the wave field, the phase singularities are points. The phase singularities show up 

>^ . at the intersection of equiphases. While turning around a singularity the phase 

I changes by 27Tp with p g Z; g = ±p is called the topological charge associated with 

■ the phase singularity and its sign is determined by the sign of the phase vortex. It 
^ ! is known that for Gaussian statistics of the field, large topological charges have a 

I small probability so that we can restrict to q = ±1 [1,2]. 

o : 

. The total topological charge present on a surface S is defined as the sum of 

O ■ the charges of the singularities located on S: Q = J^neS^^- According to Stokes' 

>► . theorem, the total topological charge Q is also equal to the accumulated phase 

I along the contour C of surface S : 

• Q^l.£ds-v^. (1) 

In this paper we study the statistics of Q and their dependence on the degree of 
disorder in a wave guide. 

In a random Gaussian speckle pattern generated by a 3D infinite medium, the 
density of singularities on a flat surface is (i = 27r/3A^ [3]. Since near field speckle 
spots are typically in size, each speckle spot contains approximatively two 
singularities. As we average over the disorder, < Q{S) >= but the fluctuations 
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of the total charge, quantified by the variance < Q'^{S) >, may be sensitive to 
both the surface S and the mean free path £ of the waves. 

Let's consider the fluctuations < Q^{R) > of the total topological charge 
contained in a circular surface S = tvR^ of radius R. One basic feature of 
< Q^{R) > is already known. If we would assume all nodal points n to have 
random charges = ±1, with equal probability and independent to each other, 
we would find that < Q^{R) >=< (^^=1 >^ ^ = TvR^d, i.e. the fluctuations 
are proportional to the surface. However, this scenario turns out to be invalid, 
at least for infinite media. Indeed it has been shown that zeros with positive 
charge tend to be surrounded by zeros with negative charge and vice versa [4-6]. 
Topological charges are not independent but tend to be screened, making the 
fluctuations grow slower than a R^ quadradic law. This screening is similar to the 
one of electrical charge in ionic fluids and plasmas. 

For a random Gaussian superposition of plane waves in space, Wilkinson and 
Freund [7] report a hnear, "diffuse" asymptotic form: < Q^{R) >^ R for large 
R. Berry and Dennis [3] use Gaussian-smoothed boundaries and show that such 
smoothed fluctuations are independent of the number of dislocations N and hence 
independent of R: < > smoothed^ cl + 0{N~^). These two methods treat the 
medium beyond distance R differently. 

In previous work, we have shown that 2D and 3D infinite media [8] reveal a 
diffuse behaviour as found by Wilkinson and Freund [7]. The role of the mean 
free path £ was also studied. In 3D, < Q^{R) > at large R depends very weakly 
on the mean free path £, with a finite value < Q^{R) > /R ioi I — oo. For 2D 
random media < Q^{R) > was seen to depend logarithmically on k£. In this paper 
we now consider a cylindrical wave guide which is a configuration used in several 
microwaves experiments by Sebbah et al [9, 10]. One additional reason for us 
to consider a wave guide is to have a genuine boundary. As a result the field is 
confined inside the wave guide so that we are sure not to forget any contribution for 
the screening process. Finally we will be able to study the influence of boundaries 
on the screening of topological charge. We study the dependence of < > with 
radius R and with mean free path £. 



2. Direct calculations 

We consider a complex random scalar wave field described by circular Gaussian 
statistics (known to be a good approximation for multiply scattered waves) [11,12]. 
The medium is a hollow conductive cylindrical wave guide of radius R and with 
infinite length, containing disorder. The problem is formulated in cylindrical 
coordinates (p, 6, z). We impose that the field derivative cancels at the boundaries. 
In the following, we first present our calculation method of the topological charge 
variance < >. 

The modes "i^rnnk (Fig. H]) of a homogeneous empty conducting cylindrical wave 
guide of radius R and infinite length are given by: 

*mnfc(p, 0, z) oc Jm ( ] exp(im6>) exp(iA:2:) with J'^{amn) (2) 
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where A: G R is the wave vector, m G Z, n G IKI and o:^^ is the n root of the first 
derivative of the m^^ Bessel function. The dispersion relation for the mode "i^mnk 
reads cj^ = k^c^ + a^^c? /E? and has a cut-off frequency ujrnn — c^mn 

c/R. The 

modes ^mnk have a mean free path rmn{k). 




^ ^ i 
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# I f f ^- 
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Figure 1. Examples of eigen modes ^mn in a cylindrical wave guide. The quantum number m 
caracterises the dependence on the angular coordinate, and the quantum number n caracterises the 
dependance on the radial coordinate: has m radial node lines and n circular node lines. 



In the presence of scattering, the averaged Green function can be written as [13]: 



f-^ J ^mn{kY - CJ^ zb le + lU/Tmn 



where ujrnn{k) is the dispersion relation of "i^rnnk Tmnk — ^mnk/c is the life time 
of the mode "ifmnk- For simplicity we shall assume that all modes have the same 
life time r — ijc although one can easily consider a more realistic model such that 
i is different for all modes. The field correlation function is proportional to the 
imaginary part of the Green function [14]. Hence we can calculate the correlation 
function between two points separated in angle by = 9' — 9 and located at same 
p and z from: 



A6>,cj) =< ^{p,9,z)^{p,9' ,z) >(x lmG{p,9,z',p,9' ,z]u) 



oc^lcos(mA0)-2^ 



T 



R 



(4) 



where the coefficients A^^(cj) are given by: 
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{{u;R/cy-al,J^+{u;RyTrr.r.r c ^ ^ 



mn 



-^"^ if ^ < amn 



(5) 

Note that for a finite mean free patli ^, the frequencies below the cut-off ujmn 
contribute as well. For sufficiently large mean free path (in FigI3] we will discuss 
how large i must be), Eq. [5] simplifies to: 




For circular Gaussian statistics, the phase derivative correlation function C^'{A9) 
can be calculated from the field correlation function using the relation [8]: 

The variance < Q'^{p) > of the total topological charge enclosed by a surface 
s{p) = Tvp^ (p E [0 : R\) with a circular contour C{p) centered inside the wave guide 
is calculated from Stokes' theorem (Eq. [1]): 

{Q\P)) = jj.P) jf ds ■ {V$(s) V$(s')) • ds' 

- A iip) r p'<^(^d^' i-.'-^'-^) - - rdAec,ip,Ae). 

C^Tr^ Jc Jc(p) de dd' I 27rJo ' 

(8) 

Eqs. [HI [71 and [IHMS] are calculated numerically. 

The results presented in Fig. [2] show that < Q^{p ^ R) > for the entire 
cross-section of the wave guide (red curve) tends to rise linearly with ujR/c. This 
global behaviour is consistent with the calculation for a 3D infinite medium (green 
curve), with a diffuse behaviour revealing charge screening. On the other hand, 
near the cut-off frequencies, < Q^{R) > exhibits sharp resonances reaching a 
maximal value of < Q^{R) >— m? for £ sufficiently large (i.e. when the limit of 
Eq. [6] is valid, figure 2a). For large m this makes < > much bigger than the 
prediction for a 3D infinite medium. We can explain the asymptotic peak value m? 
of < Q^{R) > with a simple argument. As we can see in Eq. [6l when uoR/c ~ o;^^ 
the averaged Green function is dominated by the two eigen modes ^±mn so that 
the field correlation function becomes C{9^p) ^ cos{m9) which exhibits long 
range order. Then a short analytical calculation using equations [71 and [8] shows 
that < >— m^. Note that we study the limit ^ ^ oo only with a pedagog- 
ical purpose; a truly homogneous wave guide would not exhibit Gaussian statistics. 

For a finite mean free path £ (Fig. 2b) the resonant peaks of < > are strongly 
attenuated though still visible until £ > R. The disorder introduces a coupling 
between the modes and the relative importance of the eigen modes diminishes so 
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Figure 2. (color online) Evolution of the topological charge variance < Q^{R) > enclosed by a surface 
*S = ttR^ against the product ujR/c, in red for a entire cross-section of the cylindrical conducting wave 
guide, in green dashed line for a 3D infinite medium. At the cut-off frequencies of the wave guide, sharp 
resonances arise, (a) In the limit of large calculated from Eq. [6] Figure (b) shows the damping of the 
peaks with smaller i. (b) and the zoomed inset in (a) show that the global behaviour is close to the one 

predicted for a 3D infinite medium. 



that the other modes start to contibute to the field correlation function C^. Fig. 
[3] shows the slow dependence of the peak value of the topological resonance on 
the mean free path £. Upon varying 10 orders of magnitude over the mean free 
path the fluctuations vary only by 2 orders of magnitude. The maximum value 
< Q^{R) >— m? is reached only for very large mean free path and for £ 2^ lOOi? 
the peak value is already suppressed by a factor of 10. 
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Figure 3. (color online) Semilogarithmic plot of < Q^(R) > /m^ at different dominant cut-off 
frequencies (for n = and for different m) as a function of the mean free path normalised by the wave 
guide radius. For i < R, < > approaches the value predicted for a 3D infinite mediu. Only for i as 
large as 10^, the maximum value < >~ m'^ is reached. 



In Fig. m we study how < Q'^^ip) ^ cut-off frequency Umn varies with the 

surface s = Tvp^ (with p G [0 : R]) centered inside the wave guide. This calculation 
reveals that the large value for < Q^{R) > is essentially an effect localized near 
the edge. The charges are screened in the center of the wave guide but not near 
the edges. Note that when m increases and n decreases , the unscreened charges 
localise closer to the edges. 
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(a) n = (b) m=19 

Figure 4. (color online) Evolution of the variance of the topological charge < Q'^^ip) > enclosed by a 
disk of radius p centered inside the wave guide at different cut-off frequencies ujmn and for i sufficiently 
large (a) at fixed n = and for different m (b) at fixed m = 19 and for different n. We can see that there 
is a critical radius pc- For p < pc, charges are screened and for p > pc charges of the same sign 
accumulate to give huge < Q'^nip) Pc tends to increase with m and decrease with n. 

3. Simulation 

The direct calculation of topological charge fluctuations reveals a linear depen- 
dence of the charge variance on the surface raduis caused by charge screening, 
and special frequencies for which screening seems absent. For the dominant peaks 
(n = 0), ajno ^ m and hence < Q^{R) >— m? ^ a'^^ ^ {umR/cY oc scales 
with the surface as if all the charges were independent. It is also seen in Fig. [H 
that these fluctuations mainly come from the edges of the wave guide and that the 
peak value of the resonances depends on the number of radial node lines m and 
not on the number of angular node lines n. We get deeper insight into this phe- 
nomenon using a computer simulation that generates a random circular Gaussian 
complex wave field in a transverse cross-section of the wave guide. To this end 
we write: ^{p,9) = Y^mn^rnn^mnip^O). Here ^mn{p,0) are the normalised eigen 
modes and are randomly generated coefficients obeying a Gaussian circular 
statistics with a variance fixed by: < >= Ar^n{ijj) (to be consistent with Eq.[4j). 

Fig. 5a shows the expected screening away from the cut-off frequencies. Fig. 5b 
shows that close to a cut-off frequency, charges tend to be screened in the center 
but exhibit the same sign near the edges consistently with the result obtained in 
FigH 

The twin principle [5] imposes that if the field is continuous each topological 
charge is necessarily connected by equiphase lines to a singularity of opposite sign. 
This generates the usual charge screening observed in the center, similar to infinite 
media. However in the presence of boundaries the twin principle no longer holds and 
isolated singularities may be created or destroyed at a boundary. This perturbs the 
balance between positive and negative charges. Consequently, a huge topological 
charge can appear in a circular area without an important change in the number 
of singularities. Away from the cut-off frequencies, these isolated singularities are 
independent of each other and their number is proportional to the perimeter so 
that we find < Q^{R) >— N Boundary ^ R- Howcvcr at the cut-off frequencies 
(jj = ajnnc/R the field is dominated by the weight of the two eigen modes ^^±rn,n 
with m radial nodal lines that exhibit m isolated singularities. The larger m, the 
more equiphase lines end up at the boundaries and the more singularities have 
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(c) 



Figure 5. (color online) Representation of one realisation of the random field in the wave guide, for 
different i and for different uoR/c. (•) denotes a phase singularitiy with negative charge, (o) denotes a 
phase singularity with positive charge. The colored lines denote equiphases. (a) Non cut-off frequency 
ujR/c = 20.72: Q = -2; (b) Cut-off frequency uR/c > ai9,o = 21.18227 with £ ^ oo: Q = 19] and (c) 
Same frequency ooR/c > 0:19,0 = 21.18227 but now with a mean free path £ = lOR: Q = 5 

become isolated. This increases the probabihty to have many singularities of the 
same sign, and this probability thus increases the variance < Q'^{r) >. As we can 
see in Fig. 5b these boundary singularities are not independent (which would give 
< Q'^ >^ m) but tend to be of the same sign so that < >^ m? for a large 
mean free path. The huge fluctuations at special frequencies are thus not due to 
independent charges in the total surface S but to collective effects near the edges. 

4. Conclusion 

A superposition of waves scattered by a disordered medium gives rise to a speckle 
pattern which presents a complicated network of phase vortices. We have studied 
the role of mean free path and boundaries on the screening between the topological 
charges of the phase vortices. The same linear diffuse law that relates fluctuations 
of topological charge and enclosed surface, as already found for 3i? infinite media, 
is seen. However, at the cut-off frequencies of the wave guide, giant fluctuations of 
the topological charge occur. These fluctuations are very sensitive to the disorder 
and probe the mean free path even when it is much larger than the wave guide 
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